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In this paper we investigate the effective scattering length a of spin-polarized Ne*. Due to its 

anisotropic electrostatic interaction, its scattering length is determined by five interaction potentials 

instead of one, even in the spin-polarized case, a unique property among the Bose condensed species 

and candidates. Because the interaction potentials of Ne* are not known accurately enough to 

predict the value of the scattering length, we investigate the behavior of a as a function of the five 

►7- ■ phase integrals "I>n corresponding to the five interaction potentials. We find that the scattering 

length has five resonances instead of only one and cannot be described by a simple gas-kinetic 

approach or the DIS approximation. However, the probability for finding a positive or large value of 

the scattering length is not enhanced compared to the single potential case. We find that the induced 

dipole-dipole interaction enables strong coupling between the different | J^lPMp) states, resulting in 

■ an inhomogeneous shift of the resonance positions and widths in the quantum mechanical calculation 

as compared to the DIS approach. The dependence of the resonance positions and widths on the 

input potentials turns out to be rather straightforward. The existence of two bosonic isotopes of Ne* 

enables us to choose the isotope with the most favorable scattering length for efficient evaporative 

cooling towards the Bose- Einstein Condensation transition, greatly enhancing the feasibility to reach 

, this transition. 
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Bose-Einstein Condensation (BEC) has been observed in cold dilute samples of ground state alkali atoms |l],|],|],[| 
Oh, and atomic hydrogen [pip). In 2001, the first condensate of atoms in an electronically excited state was obtained 
for metastable He[(ls)(2s) 3 Si] [f|J|], referred to as He*. All these systems have an electron configuration with only 
^vq I s-electrons in their open shells in common, resulting in an isotropic electrostatic interaction. 

. The other candidate for achieving BEC with atoms in an electronically excited state is metastable Ne[(2p) 5 (3s) 3 P2], 
referred to as Ne* in this paper. Two groups are pursuing this goal: the group of Ertmer in Hannover || and our 
in : group JnJ. Metastable neon is unique among these species in that its binary electrostatic interaction is anisotropic, 
due to its (2P)- 1 core hole Q. 

Crucial in reaching the BEC phase transition is a large ratio of "good" to "bad" collisions, i.e., a large value of the 
• elastic collision rate characterized by the total cross section a = 87ra 2 for elastic collisions with a the s-wave scattering 
length, and a small rate for inelastic collisions and other loss processes. In addition, the creation of a stable BEC 
requires a positive value of the scattering length. For metastable rare gas atoms, such as He* and Ne*, the major 
loss process is Penning ionization in binary collisions. Fortunately, the latter process is suppressed in a sample of 
spin-polarized atoms For He*, the suppression is very efficient: only spin flips due to magnetic interactions result 
in some residual ionization. Theoretical predictions and recent experimental data on residual ionization are in good 



agreement |J|bJlJ,|]i 

£Li For Ne*, the anisotropy in the electrostatic interaction determines the magnitude of the residual ionization in a 
spin-polarized gas and has a profound influence on the value of the scattering length. Theoretical estimates of the rate 
• *h , constant (3 po1 for residual ionization of Ne* predict a suppression of ionization by a factor in the range of 10 — 1000 
PlP^f , depending on the details of the interaction potentials. In experiments in Hannover and Eindhoven, a lower 
limit on the suppression of ionization by a factor of 10 has been confirmed, but so far no conclusive experimental data 
on the residual ionization rate of Ne* are available. 

In addition, the anisotropy in the interaction results in different interaction potentials Vq for the molecular states 
I J, SI) of the colliding Ne* atoms, with SI the absolute value of the projection of the total electronic angular momentum 
J = ji + 21 of the two colliding atoms on the internuclear axis. For binary collisions of spin-polarized Ne*, we have 
J=4 and SI = through 4, depending on the relative orientation of the atoms during the collision. This is illustrated 
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in Fig|l], which shows two colliding atoms in the fl — and £1 = 4 state, respectively, with the electronic angular 
momentum j'1.2 and the (2p) _1 orbital of the core hole indicated schematically. 

Because the scattering length is determined by the phase integral of the interaction potential, these different 
potentials Vq for Ne* correspond to different scattering lengths ajj. Since there is no preference for a certain relative 
orientation of the atoms (or O state) during the collision, even for spin-polarized Ne* , the elastic collision cross section 
will be determined by an effective overall scattering length a, incorporating the behavior of all five f2 states involved. 
Among the species where Bose-Einstein condensation has been achieved, the BEC candidate Ne* thus has a unique 
property. 

In this paper we investigate the relation between the effective overall scattering length a and the phase integrals 
$n of the potentials Vq. Although the potentials of Ne* are unknown at the level of accuracy needed to predict the 
value of the scattering length, it is useful to investigate the behavior of the effective scattering length as a function of 
the average value of the phase integral. In the following, we will refer to the effective overall scattering length simply 
as scattering length. 

A better understanding of the complex scattering length of Ne* is crucial in determining the feasibility of achieving 
BEC with Ne*. Important questions that need to be answered to determine the feasibility for achieving BEC are: 
(1) Is there a larger probability of encountering positive values of the scattering length, as compared to the 75% 
probability for the single potential case? (2) What is the probability for finding a sufficiently large elastic total cross 
section for efficient evaporative cooling? (3) How does the availability of two bosonic isotopes of Ne*, 20 Ne* and 22 Ne* 
(with a natural abundance of 90 % and 10 %, respectively) influence these chances? 

The values of the rate constant for residual ionization f3 pcl and the scattering length a which we need to achieve BEC 
with Ne* for typical experimental conditions in our experiment, are summarized in Fig. |j. This so-called feasibility 
plot for achieving BEC with Ne* shows the calculated number of atoms N c with which quantum degeneracy is achieved 
as a function of (3 po1 and the absolute value of the scattering length |o|. The calculation of this feasibility plot is based 
on the kinetic model of Luiten et al |i~7|], including trap losses. The criterion for crossing the BEC transition has been 
set to iV c > 1 x 10 5 atoms. The initial conditions of the evaporative cooling process are taken equal to experimental 
values which we are able to produce, i.e., N — 1.5 x 10 9 , T = 1.2 mK, r — 8 s and r\ = 5.5, with N the number of 
atoms, T the temperature, r the lifetime of the atom cloud, 77 = E/iksT) the truncation parameter, E the energy 
and ks the Boltzmann constant. 

From Fig. || it is clear that BEC of Ne* is feasible for (3 P ° l < 2.5 x 10 -13 cm 3 s _1 (a suppression of 500) and a c > 75 
ao, with a the Bohr radius. Therefore, we calculate the probability P c that either of the bosonic isotopes of Ne* has 
a total cross section larger than a c = 87ra 2 = 1.4 x 10 5 a 2 ,. A larger suppression of ionization of course allows for a 
smaller value of the scattering length and vice versa. 

As a first order estimate, in a simple gas-kinetic picture, one might expect the effective elastic total cross section 
cr, referred to as elastic cross section in the remainder of this paper, to be a weighted average over the elastic cross 
sections an — 8™^ of the different molecular states. Applying this approach to the anisotropic Ne* problem results in 
a large enhancement of the probability of encountering large values of a as compared to all systems with an isotropic 
interaction potential. Moreover, the elastic cross section is always larger than a rather large lower limit. 

A more sophisticated approach can be found in the Degenerate Internal State (DIS) method |2^]. In this method, 
the energy splitting of the internal states is neglected and a is given by the weighted average of the contributing 
scattering lengths clq. For cold collisions of hydrogen atoms, the DIS method results in values for the scattering 
length that compare well to the outcome of full quantum calculations |lg|]l9|| . 

The definitive approach to determine the scattering length a of Ne*, of course, is a full five-channel quantum- 
mechanical calculation. Although this numerical approach supplies the correct answer to our problem, it has the 
disadvantage that the results are not always easy to understand in terms of the properties of the input potentials. We 
use the results of the numerical calculation to check the validity of the different analytical approximations described 
above, which in general give more insight. 

This paper is organized as follows: first, the available ab-initio potentials and the calculation of t he diff erent 
phase integrals <&n are discussed in Sec. [Il[ In Sec. Ill the s ingle potential scattering length (Sec. Ill A ), the 



gas-kinetic approach to d efining an elastic cross section a (Sec. IIIB) and the scattering length obtained with the 
DIS approximation (Sec. QIC) are discussed. We then present the results of our quantum mechanical, numerical 
scattering calculation in Sec. [V. To end we present our conclusions in Sec. [v| All calculations are performed for 



both bosonic isotopes of Ne*, ^Ne* and 22 Ne 
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II. INTERACTION POTENTIALS 



There are five adiabatic molecular states | J, O) that connect to the spin-polarized Ne* + Ne* asymptotic limit with 
total electronic angular momentum J = 4 and total spin 5 = 2. The degeneracy of the O = state is f , that of all 
others is equal to 2. As input potentials for Ne* we have used the short-range ab-initio potentials of Kotochigova et al. 
|0|,fnj, which are available in the range R < Ri, with i?i = 60 ao for = 4 and R\ = 120 ao for all other potentials. 
Typical values of the well depth e and its position R m are e ~ 30 meV and R m « 10 ao- The long-range behavior of the 
potential curves is dominated by the attractive induced dipole-dipole interaction —Cq/R 6 . The ab-initio potentials 
have within < 3 % identical Ce coefficients, since the long-range interaction is dominated by the (3s) valence electron 
0,0. We use a single Cq coefficient with a value of Cq — 1938 a.u. as calculated by Derevianko et al fl(|, based on 
the static polarizability of neon. 

We characterize the potentials Vq by their classical phase integral 



I* Its noo 

®n= k n (R)dR+ / k n (R)dR 

= (1) 



with kn(R) the local wave number and R c the classical inner turning point for zero collision energy. We choose R s 
such that for R < R s the energy splitting AVn t n'(R) = Vn(R) — Vq'(R) between the O-potentials dominates over the 
rotational coupling AV ro t{R) = -[ti 2 /(2[iR 2 )]8\n-n>\=i V( P ( P + = V( J ( J with A 4 the reduced 

mass, I the rotational angular momentum and P = J + I the total angular momentum. For R > R s the opposite 
holds. 

The first part $>q <Rs of the phase integral has been calculated by numerical integration in the interval [R c , R s ], with 
R s = 20 ao. The contribution for R > R s to the phase integral is calculated analytically assuming a pure long-range 
-C e /R 6 behavior, resulting in $ R>R ° = 3.34 ir. 

The numerical results are given in Table | for both bosonic isotopes 20 Ne* and 22 Ne* . We see that the different 
O-states have very different values of $n, varying by as much as A$o.si' = &n — = 0.8 tt. This implies both 
a different number of bound states and different positions of the resonances in ao- Asymptotic behavior of, or a 
resonance in the scattering length occurs when a quasi-bound state lies close to the dissociation limit or has just 
moved into the continuum. 

Because the interaction potentials of neon are not known accurately enough to predict a, we have to vary $q over 
a range equal to tt to predict the range of a values that we can expect for the spin-polarized Ne* system. Therefore, 
we introduce a scanning parameter 4> that we add to $n to create a modified phase integral 0q according to 

<t>n = $n + <t>> (2) 

with 4> G [0,7r]. 

For simplicity, we assume for now that the phase differences between the O potentials A$si,f2' are constant and 
equal to the ab-initio values given in Table Q. Later on in Sec. IV, we will also vary A$f2,n' over an interval tt to 
investigate the influence of A$n,n' on the scattering length a. 

The classical phase integrals of the two bosonic isotopes of Ne* are related by the mass scaling rule 

22 fe = (22/20) 1 / 2 2 Vn- (3) 

Using the average phase integral over all O states ( < i ) )o = 16.5 tt, we find an isotope shift equal to 0.81 tt. This simple 
relation between the phase integrals of 20 Ne* and 22 Ne* enables us to compare very easily the results obtained for 
20 Ne* with those for 22 Ne*. 



III. ANALYTICAL APPROACH 

A. Single potential scattering length 

The semiclassical analysis of the scattering length in atomic collisions by Gribakin and Flambaum |2^] yields for 
the s-wave scattering length of a potential with a long-range behavior —Cq/R 6 

7T 

ao = a bg [l - tan(0n - ^)] 
o 



a bg = cos(^/4) [ v ^Z^/4n] 1 / 2 [r(3/4)/r(5/4)], (4) 
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with r() the Gamma function and cib g the background value of the scattering length. The latter is fully determined 
by the long-range behavior of the potential and is equal to a,b g = 44.3 ao- The position of the resonance in is equal 
to c/)q s mod it = it/2 + tt/8 = 5ir/8. We define the width of a resonance in an as 

Tj2 = ^n(3a feg ) - 4>n(-ab g ), (5) 

around the resonance position <p r ^ s . The scattering length varies around ab g with a probability of 75 % of being 
positive. The probability that at least one of the bosonic isotopes of Ne* is positive is much larger: 94 %. 

B. Gas-kinetic model 

In a simple gas-kinetic approach, we define the elastic cross section as a weighted average of the elastic cross sections 
0-fi 

a = (8™^) 
J 

= a bg Y, w n [l-tan(0n-^)] 2 . (6) 
o=o 

Here dbg — 87ra^ = 0.5 x 10 5 a 2 , is the background value of the elastic cross section and wq, is the amplitude of the 
projection of the initial asymptotic rotational state on the basis, with u>o = 1/3 and u>n>o = \/2/3. We assume 
that we are in the low temperature limit where the condition h«l holds, with k — ^/2[iE jfi the asymptotic wave 
number and E the collision energy in the reduced system. 

In Fig. H we show the results for the elastic cross section of 20 Ne* (solid lines) and 22 Ne* (broken lines) as a function 
of <fi € [0, 7r] (Eq. (0)). Two important characteristics in the elastic cross section of Ne* immediately catch the eye. 
First, we see a rather large minimum value <J m in for the elastic cross section and an increase of almost a factor of two 
in the probability P c that the elastic cross section is large enough to make BEC of Ne* feasible as compared to the 
single potential case, as can be seen in Table [n]. 

Secondly, we see five resonances in a, attributable to the five different O-states of Ne*. This characteristic behavior 
does not depend very much on the specific values of A<I>si,o as long as they are not very small (< 0.05 it). This 
behavior is very different from the single potential case, where we can encounter an elastic cross section equal to zero 
and there exists only one resonance. 

The general picture is the same for both isotopes. However, the elastic cross section of 22 Ne* is shifted with respect 
to the elastic cross section of 20 Ne* by the isotope shift of 0.81 it. Depending on the actual phase integral <t>n of the 
20 Ne* system, it can thus be advantageous to switch to the less abundant bosonic isotope 22 Ne*, to optimize the value 
of the elastic cross section. Choosing for each value of <p the isotope with the largest a yields an even larger minimum 
value of a, as can be seen in Table [n} In addition, P c increases from 78 % in the single isotope case to 99 % for either 
isotope. 

C. DIS model 

Next, we investigate the scattering length in the DIS approximation, which has proven to be quite insightful for 
hydrogen [|lg| , |l9|| . In this approach, the energy splitting of the internal states is neglected. For Ne*, this means that 
the rotational splitting AV ro t between the partial waves is neglected. In the DIS approximation, the scattering length 
is given by a weighted average of the five asi's involved 

a = (an) 

4 

= a bg ^ ~ ta.n.((f>n ~ -)}. (7) 
n=o 

From Eq. (Q) it is clear that the resonance positions (/>ff s of a in the DIS approximation, coincide with the single 
potential resonance positions. They are completely determined by the values of A<&n.o 



,5 

y f2 ' ' 



t-res _ _ A $n i0 ) mod 7T. (8) 
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Figure || shows the scattering length a as a function of the scanning parameter <fi £ [0, w]. Again we observe five 
resonances in the scattering length. In this figure we have also plotted the behavior of the single potential scattering 
length dn=4 (broken line), showing clearly that the single potential resonance positions coincide with the resonance 
positions of a in the DIS approach. 

Taking the weighted average of an does n ot cha nge the total probability for a positive scattering length (75 %) as 



compared to the single potential case (Sec. [Ill A| ) . However, the probability of encountering a large value of a does 
increase, as can be seen in Table |n| Both the width of the resonances (Eq. <^)) and the derivative da/ dip of a with 
respect to cf> determine this probability. From Fig. |] it is clear that the total width T of all five resonances combined 
is much larger than the single potential resonance width. The width of each resonance is not only determined by the 
relative weight wn of its single potential scattering length but also by its relative position <jf^ s in respect to the other 
resonances, and therefore depends sensitively on the phase differences A^n^' between the potentials. An increased 
width and derivative of a in the DIS approach as compared to the single potential case therefore lead to a much larger 
P c (Table [n]) . Choosing the most advantageous bosonic isotope improves again a m i n and P c as well as the probability 
for encountering a positive value of a (95 %) as compared to the single isotope case (Table H). 

IV. QUANTUM MECHANICAL NUMERICAL CALCULATION 

A. Model 

For the Ne* system, we can distinguish two regions of interest in the potentials Vq(R) (Fig. ||). At small internuclear 
separations (region I), i.e. R < R s , the splitting AVq^'(R) between the interaction potentials Vq is larger than the 
differential rotational coupling AV ro t(R)- As a result, there is no coupling between different f2 states but coupling 
between different / states. In this region 17 is a good quantum number and the \JttPMp) basis is the proper 
representation. Here Mp is the projection of P on the quantization axis. Both P and Mp are conserved during the 
collision. In a semiclassical picture we can describe this regime in the following way. The large fi-splitting results in 
a rapid precession of J about the internuclear axis, much faster than the rotation of the axis itself as determined by 
the value of \l\. The projection of J on the internuclear axis is thus conserved, while the magnitude of / changes due 
the changing orientation of J with respect to the space-fixed total angular momentum vector P. 

At large internuclear separations (region II in Fig. |^), i.e. R > R s , the rotational coupling dominates the interaction: 
the splitting AVh,n'(i?) is smaller than the differential rotational coupling V rot (R). The relative motion of the colliding 

atoms results in a rotation of the internuclear axis with respect to J and therefore in a change in the value of fl (Fig|l|). 
In this region fl is not a good quantum number but I is, and we use the | JlPMp) representation. 

Ultra-cold collisions between spin-polarized Ne* atoms are described by a five-channel problem: five tt= 0..4 
channels in region I and five 1 — 0,2, ..8 channels in region II. Only even partial waves contribute due to the symmetry 
requirement of the wave function for bosons. The rotational energy barrier (5.6 mK for 1=2, located at 78 ao) is 
always much larger than the collision energy (< 0.5 mK). For this reason, higher order partial waves (/ ^ 0) do not 
contribute to the incoming channel. However, in region I, the short-range interaction AVo.fi' still couples the single 
incoming channel |J = 4Z = 0P = 4 Mp = 4) to higher order partial waves. Because the tunnelling probability 
for I is negligible (< 10 -5 ), they only contribute to the elastic scattering process when they again couple to the 
|J = 4Z = 0P = 4 M P = 4) initial state. 

In our calculations we assume that the intermediate region, where V ro t(R) and Vh,c(P) are of the same order 
of magnitude, is arbitrarily small, i.e., we assume a sudden transition from region I to region II at R = R s . The 
scattering problem now reduces to potential scattering and we can solve the uncoupled problem in region I in the 
| JOPMp)-basis and in region II in the | J7PMp)-basis, for each channel. After transformation of the solution ui(R) 
in region I at R = R s from the | JflPMp)-basis to the | J7PMp)-basis, we connect it continuously to the long-range 
solution uji(R) at R s , assuming equal local wave numbers pS" 



In very good approximation |24 25J, we can summarize the behavior of the atoms in region I by means of the 



accumulated phase method. In R — R s the radial wave function of a single uncoupled channel is then given by 

4 2 ( J R s )=sin(^<^+0 + 7 r/4), (9) 

with the extra phase shift 7r/4 due to our choice of using the classical phase rather than the quantum mechanical 
accumulated phase. After transformation to the \ JlPMp) basis the 5x5 solution matrix is given by 

u\{R s ) = T lQUl (R s ), (10) 

with Tin the elements of the transformation matrix T between the \JVtPMp) and the \ JlPMp) basis. 
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In region II, the evolution of the radial wave function u l H for R > R s is governed by the radial Schrodinger equation 



9R2 ^(R)+^-^-^)u> II (R) = o. (ii) 

The connection of the short (I) and long-range (II) solutions in the \ JlPMp) representation at R — R s determines 
the boundary conditions for the numerical integration of u\ I from R s to R a , where the asymptotic limit of a vanishing 
potential is valid. At R a , the numerical solution is connected to the asymptotic radial wave function 

U l n (R a ) = [ Aie -i(kR a -l*/2) + Bie i(kR.-ln/2)^ ^ 

to determine the scattering matrix S = B A~ ■ From the scattering matrix we then obtain the scattering length a 

a= _ lim tan(ln(5oo)/2^ 
fc^o k 



B. Results 

Using the ab-initio potentials of section p we h ave calculated the scattering length of 20 Ne* using the quantum 
mechanical calculation described in section [V A| . Again, we have performed these calculations for <f> 6 [0, it], to 
determine the range of scattering length and elastic cross section values that we can expect for the Ne* system. 
Figure H shows the scattering length a (solid line) as a func tion of tfi. For comparison, we have included the result for 



a obtained with the DIS model (broken line) (Subsec. Ill C ), in which the positions of the resonances in a correspond 



to the single potential resonance positions in ciq. The resonance positions 0^ es of a are shifted in an inhomogeneous 
way with respect to those obtained with the DIS method and also the widths Tq of the individual resonances differ 
significantly. We can no longer attach an fi-label to each of the resonances, as is possible in the DIS approach. In 
the DIS model, we neglect the rotational splitting V ro t between the partial waves, so that effectively all partial waves 
are equivalent. Apparently, both the resonance positions as well as the widths of the individual resonances in the 
numerical results for a are influenced by the coupling to higher order partial waves. This is not surprising: close lying 
bound states from other ^-potentials will most likely shift the bound state and thus the corresponding resonance 
position. Clearly, neglecting the rotational splitting of the internal states is not justified in the case of Ne*. 

However, despite these qualitative differences in the behavior of a, the quantitative behavior is the same as in the 
single potential case. The probabilities for encountering a positive a value and an elastic cross section a c are equal 
to those found in the single potential case (Table ||). Although the total width of the resonances is larger, this is 
compensated by a decrease in the derivative of a. Moreover, choosing the most advantageous bosonic isotope for each 
value of the scanning parameter 4> leads to a similar increase in these values for both the full quantum mechanical 
calculation as the single potential case. 

To investigate the influence of A$n,a' on the positions of the resonances in a, we have varied the classical phase 
difference A$4 ! o of one of the £7 potentials (£1 = 4), while keeping the other classical phase differences fixed at their 
ab-initio value (Table |). Starting at its ab-initio value, the classical phase difference 

A«£ 4 ,o = A$4,o + A0, (14) 

with A(f> £ [0,7r], is varied over A</>4.o G [—0.54 7r,0.46 tt]. In this way, the bound states in the £1 = 4 potential 
encounter the bound states in all other £1 potentials. The position of the resonances <pQ B for this modified set of 
ab-initio potentials is determined in the usual way, by scanning the parameter <fi in the phase integral (f>n (Eq. (||)) 
over the range [0, n]. 

The results are shown in Fig. ^, where we have plotted the position <f>Q S of the five resonances in a as a function of 
the shift A<f> in the bound states of the £1 = 4 potential. The broken lines are drawn to guide the eye. We observe that 
two of the resonances remain at a fixed position, while all three others shift proportional to A</> with a slope equal 
to (-0.34 ± 0.02)tt and separated by (0.30 ± 0.02)vr. Narrow avoided crossings occur when 'constant' (j)^f s meet 4>^ s 
varying like cx —Acf> at Acj) = (0.46 + n)ir and A<fi = (0.52 + n)-7r, with n = 0, 1, 2, ... In addition, very broad avoided 
crossings occur between 0^ es varying like oc — A0 at A</> « (0.1 + n)7r, with n = 0, 1, 2, ... We have plotted <jf^ s over a 
range A<p £ [0, 2ir), to show the avoided crossings. Varying one of the other classical phase differences A</>q^4.o yields 
similar results. 
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Clearly, the quasi-bound states of the system are no longer pure \JQPMp) states. Apparently, the quasi-bound 
states whose resonance positions (jf^ s vary proportional to Acf> are strongly coupled to the \JQ = APMp) state and 
those whose <j>Q S remain constant are only very weakly coupled to the \Jft = APMp) state. When two strongly 
coupled resonances approach each other, a strong, broad avoided crossing occurs. Similarly, a weak, narrow avoided 
crossing occurs when a weakly coupled resonance approaches one of the other resonances. 

This picture is consistent with the behavior of the width of the resonances as a function of A0. This becomes 
clear when we look at a simpler, more transparent two-channel, i.e. J = 1, numerical calculation of the scattering 
length. The results of this calculation are consistent with the full five-channel calculation. In Fig. ^a) we have plotted 
the phase integrals (jf^ at which a resonance in a occurs as a function of A</> for J = 1. Both resonance positions 
vary proportional to A</> with a slope equal to (—0.50 ± 0.02)7r, and with broad avoided crossing between them at 
A<p = (0.1 + n)ir, with n = 0, 1, ... They are strongly coupled, which is also reflected in the behavior of the width of 
the resonances Tq (Figure ||c)), which varies as a sine (broken line) between and T. At the avoided crossings the 
widths of the resonances become equal. The total width of both resonances combined r = Tq + Ti is conserved, but 
one resonance is wide while the other is narrow. This periodic change in the resonance width is due to the periodic 
change in the coupling of both quasi-bound states to the incoming I = channel with changing 0Pf s . 

The weakly coupled case can be illustrated by reducing the value of Cq, because in the absence of an induced 
dipole-dipole interaction —Cq/R 6 no coupling between the quasi-bound states is possible. This can be understood in 
the following way. With a decreasing dipole-dipole interaction the rotational energy barrier increases both in height 
and width. As a result, the region around R s in which rotational coupling between the different £1 states takes place 
decreases and eventually vanishes. Figure || b) and d) show the resonance positions and widths for a vanishing value 
of Cg, simulated by assuming a modified value C' & — Cq/100. One resonance position remains constant, while the 
other varies proportional to A<p and the widths of the resonances are small and remain constant, except at the avoided 
crossing where they approach each other. This is the behavior seen for some of the resonances in a for J = A and 
is qualitatively the same as the behavior of the resonance positions and widths in the DIS model, where coupling 
between the quasi-bound states is not taken into account. 



V. CONCLUDING REMARKS 



Elastic collisions between spin-polarized Ne* atoms are governed by multiple interaction potentials. This unique 
property of Ne* among the BEC species and candidates is a result of the anisotropic interaction between them. Both 
simple analytical and full quantum mechanical calculations of the scattering length a of Ne* show that the resulting 
scattering length has five resonances. A simple gas-kinetic picture yields very favorable but unrealistic results for the 
clastic collision cross section that are not compatible with the numerical calculations. This approach is only valid for 
an incoherent mixture of £1 states. Comparison between the numerical results and the DIS model reveals that a is 
also not simply a weighted average over the single potential resonances an and that the resonances in a cannot be 
assigned to a single ft state. Although the DIS approach assumes a coherent mixture of SI states, coupling between 
the quasi-bound states is not taken into account, and it therefore does not describe the Ne* system accurately. The 
overall behavior of a is similar to that of the usual single potential scattering length: neither the probability for 
encountering a positive or a large value of a is enhanced by the presence of five instead of one resonances. 

The presence of an induced dipole-dipole interaction leads to strong coupling between the different J! states and 
causes a broadening of the resonances, resulting in quasi-bound states that are a linear combination of different 
| JVlPMp) states. This coupling between the different | JQPMp) states in turn leads to the inhomogeneous shift of 
the resonance positions and widths in the quantum mechanical calculation as compared to the DIS approach. However, 
the dependence of the resonance positions and widths on the input potentials is quite straightforward. The resonance 
positions vary either directly proportional to the relative phase differences between the fi potentials or not at all, 
depending on the exact composition of its quasi-bound state. The width of the strongly coupled resonances (whose 
positions vary oc — A<f>) is determined by the coupling to the ingoing channel, which varies periodically with Atp. The 
width of the weakly coupled resonances (whose positions do not vary with Acf>) is constant, since the coupling to the 
ingoing channel does not change. The total change in all five resonance positions is always equal to — A<j> and the 
total width is conserved. 

The possibility to choose between the two bosonic isotopes of Ne* to optimize the value of the elastic cross section, 
greatly enhances the prospects for achieving BEC with Ne* (Table |lj) . Large beam fluxes of both bosonic isotopes, 
crucial in obtaining favorable initial conditions for efficient evaporative cooling, have been realized at the Eindhoven 
experiment |26|, therefore choosing the isotope with the most favorable scattering length is feasible. 
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FIG. 1. Schematic view of two colliding Ne* | 3 P2) atoms in a spin-polarized gas (S = 2, J = 4) for both the f2 = and 
the SI = 4 state. The orientation of the electronic angular momentum ji^ and the (2p) _1 core hole of the individual atoms is 
indicated schematically. 
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FIG. 2. Feasibility plot for reaching the BEC transition with Ne* in the Eindhoven experiment, showing the number of 
atoms N c with which quantum degeneracy is achieved as a function of the rate constant j3 v ° for residual ionization and the 
absolute value of the scattering length \a\. The broken lines correspond to a c and (3^° l for which BEC of Ne* becomes feasible. 



9 




FIG. 3. Elastic cross section a = (87ran) in a gas-kinetic approach for 20 Ne* (solid line) and 22 Ne* (broken line), as a 
function of the phase integral <f> £ [0, n]. We observe five resonances due to the five contributing potentials, which also results 
in a minimum value <j m in > for a. 




FIG. 4. DIS result for the scattering length a (solid line) as a function of the parameter cf> £ [0, it], showing five resonances that 
are labelled with an SI- value because they are located at the position of the single potential resonances in an. For comparison 
we show the behavior of 0,4 (broken line). 
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FIG. 5. Interaction potential Vo for the Q, — molecular state of two colliding spin-polarized Ne* atoms, with both range 
I (dominant fi-splitting) for R < R s = 20 ao and range II (dominant rotational energy splitting) for R s < R < R a — 200 ao 
indicated. 




FIG. 6. Full quantum mechanical calculation of the scattering length a (solid line) as a function of 4> £ [0, 71"]. For comparison 
we also have depicted the scattering length a obtained with the DIS model (broken line). Both position and width of the 
resonances in a differ from the DIS result: an fl-label cannot be attached to any separate resonance. 



11 



0.0- 



0.0 0.5 1.0 1.5 

A(|> (k radians) 



2.0 



FIG. 7. Numerical calculation of the phase integrals 4>n a (filled squares) at which a resonance in a occurs for a modified set 
of potentials, with the phase integral difference varying from its ab-initio value —0.54 ir to +1.57T by varying A<f> over 

2tv. In this way, the bound levels of the = 4 state 'encounter' the bound levels in all other fi-states. The broken lines are 
drawn to guide the eye. 
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FIG. 8. Numerical two-channel (J = 1) calculation of position (a and b) and width (c and d) of the resonances in a for a 
modified set of potentials, with A0i,o varying from its ab-initio value —0.61 7r to +1.39 7r by varying A<f) over 2n. The figures 
on the right hand side (b and d) show the results for a reduced C' & = Ce/100. The broken lines are drawn to guide the eye and 
the arrows indicate the avoided crossings. In the presence of an induced dipole-dipole interaction (a) the resonance positions 
vary both proportional to A(f) and are strongly coupled, and the width of the resonances (c) varies as a sine. For a small 
dipole-dipole interaction (b) one of the resonance positions remains constant, while the other varies proportional to A<f>. The 
width of the resonances (d) remains constant, except at the avoided crossing where they coincide (arrow). 
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TABLE I. Classical phase integral $n and its difference Af&j^o = $n — $0 with respect to the Q = potential of the 
spin-polarized adiabatic molecular fi-states, connecting to the Ne* + Ne* asymptotic limit with J= 4 and S= 2. The O states 
are labelled by Q g , where the gerade label g reflects the symmetry of the electron wave function under inversion around the 
center of charge. Data are given for both bosonic isotopes 20 Ne* and 22 Ne* of Ne*. 
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TABLE II. Minimum value <j m in of the elastic cross section and probability P c for a > a c — 1.4 x 10 5 Oq for a single isotope 
( 20 Ne*) and the set of two bosonic isotopes 20 Ne* and 22 Ne*, as calculated with the single potential semiclassical model, the 
gas-kinetic model, the DIS model and the quantum mechanical numerical calculation. 



model 


O'mi 


single isotope 
n (units of 10° ao) 


Pc (%) 


O'mi 


either isotope 
„ (units of 10 5 a 2 ) 


Pc (%) 


single potential 







42 




0.16 


61 


gas- kinetic 




0.7 


78 




1.3 


99 


DIS 







72 
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